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Three-jet event-shape distributions can be exploited to investigate the dependence of 
hadronisation effects on the colour and the geometry of the underlying hard event. We 
present here the first comparison of data in e + e~ annihilation and state-of-the-art the- 
oretical predictions, including resummation of large logarithms at next-to-leading loga- 
rithmic accuracy matched to exact next-to-leading order and leading non-perturbative 
power corrections. 

1 Power corrections to multi-jet event shapes 

The remarkable success of the QCD description of two-jet event shape distributions has made 
these observables one of the most useful tools to test our understating of the dynamics of 
strong interactions, both in the perturbative (PT) and non-perturbative (NP) regime. This 
is because event-shape distributions span a wide range of physical scales, from the region 
where the event shape V is large, described well by fixed-order QCD, to the exclusive V — > 
region where hadronisation effects dominate, through the intermediate region where one 
needs to resum large infrared and collincar logarithms. The combination of next-to-leading 
order (NLO) predictions and next-to-lcading logarithmic (NLL) resummation, supplemented 
with non-perturbative (NP) hadronisation corrections provided by Monte Carlo (MC) event 
generators, has lead to one of the most precise determinations of the QCD coupling a s [jQ. 

In view of the fact that hadronisation corrections are suppressed by inverse power of the 
process hard scale Q, in recent years it has been attempted to describe two-jet event shape 
distributions at hadron level by simply adding to the NLL resummation the NP shift (SV) 
originated by leading 1/Q power corrections, which is a reliable approximation as long as 
(SV) <C V. The shift has a remarkably simple structure, being the product of a calculable 
coefficient cy, which depends on the considered shape variable, and a genuine NP quantity 
(fc t ) Np , the mean transverse momentum of large-angle hadrons produced in the collision, 
which is variable independent (universal). The universality of (fct) NP , and hence of 1/Q 
power corrections, has been thoroughly tested both in e + e~ annihilation and DIS, and is 
found to hold within 20% (see [5] for a recent review). 

This universality property is based on two facts. The first is that particles responsible 
for leading power corrections are low transverse momentum hadrons in a central rapidity 
region, away from the hard jets. Any of these hadrons k contributes to a two-jet event shape 
V with an extra SV(k)Q ~ ktfv(v)i with kt and r\ the hadron transverse momentum and 
rapidity with respect to the jet axis. The second is that central hadrons are distributed 
uniformly in rapidity. This ensures that in the region (5V) <C V , where only leading power 
corrections are important, the dependence of (SV) on rapidity and transverse momentum 
gets factoriscd [3]: 



Among all models that, for two-jet events, predict a uniform rapidity distribution of 
central hadrons, the dispersive DMW approach |4] makes it possible to extend eq. (JTJ) to 
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multi-jet event shapes, where there is no natural way to identify kt and r\. The starting 
point is the probability dw(k) of emitting a soft dressed gluon k from a quark-antiquark 
pair (whose momenta are p and p) in a colour singlet: 

dw(k) = C FTT dr]- , rj = - In — , fc t = — , 2) 

where a s is the physical CMW coupling [5]. The CMW coupling is then extended at low 
transverse momenta via a dispersion relation, and the very same probability dw(k) is ex- 
ploited to compute NP corrections [I]. The resulting shift (SV) has the same form as in 
cq. ([T]), where the cy coefficient is identical and the NP parameter (fc t ) NP can be expressed 
in terms of ao(w)> the average of the dispersive coupling below the merging scale fj,j, as 
follows [6]: 

An r^i Ah 

(k t ) NP = -^-C F M (a (fii) - a s {Q) + 0(a 2 s )) , a (fn) = / — a s (k) . (3) 
n Jo M/ 

Here the Milan factor M. accounts for non-inclusiveness of event-shape variables. 

One can now naturally extend the above analysis to multi-jet event shapes, where the 
soft dressed gluon probability is given by 



X£j K ij 2tt TT 2 Pl k 2p lPj 



with Ti the colour charge of hard parton pi, and Kij and r\ij the invariant transverse mo- 
mentum and rapidity with respect to the emitting dipolc ij. This gives the following result 
for the shift: 

( 5V ) = ^T M - + °^s)) E(- f ' ' f i) c v ] ■ ( 5 ) 

i<j 

The above expression states that NP corrections to multi-jet event shapes depend on the 
same parameter ao(/i/) encountered for two-jet shapes. Moreover, they depend in a non- 
trivial way on the colour of the underlying hard event through the correlation matrices Tj • Tj 
and on the event geometry (the angles between the jets) through the calculable coefficients 



c 



0. 



(ij) 
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The simplest environment in which the validity of cq. ([5]) can be tested is three-jet events. 
Here colour conservation ensures that the colour matrices T$ ■ Tj are in fact proportional to 
the identity, thus simplifying considerably both the PT and the NP analysis. 

2 Results for three-jet event shapes in e + e" annihilation 

Two three-jet event shapes have been studied in e + e~ annihilation, the _D-parameter [S] 
and the thrust minor T m [5]. Both variables are small when the three hard jets are in a 
near-to-planar configuration, and measure radiation outside the event plane. 

We present here the first ever comparison of theoretical predictions for D and T m differ- 
ential distributions and existing data provided by the ALEPH collaboration [10] . Theoretical 
predictions are at the state-of-the-art level, that is NLL resummation matched to the NLO 
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calculation obtained with nlo JETH — h [TT] , and leading 1/Q NP corrections computed with 
the dispersive method [HUH]- Events with three separated jets are selected by requiring the 
three-jet resolution parameter 1/3 in the Durham algorithm to be larger than y cu t- It is then 
clear that different values of y C ut correspond to different event geometries. 

Figure Q] shows the result of a simultane- 
ous fit of a s (Mz) and ao(/J,i = 2GeV) for the 
D-parameter distribution at Q = AIz cor- 
responding to j/cut = 0.1 and y cut = 0.05. 
The 1-ct contour plots in the a s -ao plane 
are plotted together with results for other 
distributions of two-jet event shapes. There 
is a remarkable consistency among the var- 
ious distributions, thus strongly supporting 
the idea that universality of 1/Q power cor- 
rections holds also for three-jet variables. 
This leads to the non-trivial implication 
that leading power corrections are indeed 
sensitive to the colour and the geometry 
of the hard underlying event, and more- 
over this dependence is the one predicted 
by eq. (0). 
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Figure 1: Contour plots in the a s -ao plane 
for the _D-paramctcr differential distributions 
corresponding to two different values of y C ut- 

The comparison to data is less satisfac- 
tory for T m , as can be seen from Fig. [2j 
There one notices a discrepancy between 
theory and data at large values of T m . To 
track down the origin of the problem, one 
can look at hadronisation corrections pro- 
duced by MC programs, defined as the ra- 
tio of the MC results at hadron and par- 
ton level. From the plots in [10] one can 
see that hadronisation corrections for the D- 
parametcr are always larger than one, corre- 
sponding to a positive shift, consistent with 
our predictions. On the contrary, hadro- 
nisation corrections for T m become smaller 
than one at large T m , a feature that will 
never be predicted by a model based on a 
single dressed gluon emission from a three 
hard parton system. This issue is present 

also in the heavy-jet mass and wide-jet broadening distributions, and requires further theo- 
retical investigation. 




Figure 2: Theoretical predictions for T m dis- 
tribution plotted against ALEPH data for 
three different values of y cu t- 



3 Extension to other hard processes 

Observables that measure the out-of-event-plane radiation in three-jet events can be intro- 
duced also in other hard processes. 

In DIS two observables have been already measured. One is a variant of T m [12], where 
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all momenta are in the Breit frame, and the event plane is formed by the virtual photon 
direction and the thrust major axis, defined as the direction that maximises the projection of 
transverse momenta. Differential T m distributions have been measured both by the HI [13| 
and ZEUS collaboration [14], and fits of experimental data are currently in progress. The 
other observable is the distribution in the transverse energy correlation EtEtC(x), defined 
as [H] 



The interesting features of the EtEtC(x) distribution arc that it approaches a constant for 
small x and that it has fractional power corrections. 

In hadron-hadron collisions one can consider for instance the production of a Z boson q 
in association with a hard jet pj e t- The event plane is determined by the beam and the Z 
direction, and one can study [16] 



where the sum is extended to all hadrons not too close to the beam pipe, and the normal- 
isation is fixed so as to cancel systematic uncertainties in the jet energy scale. In order to 
compare data with existing predictions, rjo should be taken as large as is experimentally 
possible. The interest in this variable is that its distribution can take large corrections from 
the underlying event, thus making it a useful tool to tune MC models of minimum bias and 
multiple hard collisions. We look forward to experimental investigations in this direction. 
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